Spinon Phonon Interaction and Ultrasonic Attenuation in Quantum Spin Liquids 
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Several experimental candidates for quantum spin liquids have been discovered in the past few 
years which appear to support gapless fermionic S — | excitations called spinons. The spinons may 
form a Fermi sea coupled to a U(l) gauge field, and may undergo a pairing instability. We show that 
despite being charge neutral, the spinons couple to phonons in exactly the same way that electrons 
do in the long wavelength limit. Therefore we can use sound attenuation to measure the spinon 
mass and lifetime. Furthermore, transverse ultrasonic attenuation is a direct probe of the onset of 
pairing because the Meissner effect of the gauge field causes a "rapid fall" of the attenuation at T c 
in addition to the reduction due to the opening of the energy gap. This phenomenon, well known 
in clean superconductors, may reveal the existence of the U(l) gauge field. 
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Quantum spin liquid in dimensions greater than one 
is a long sought state of matter which has eluded ex- 
perimental investigation until recently. [l[ We define the 
quantum spin liquid as an insulator with an odd number 
of electrons per unit cell which does not order magnet- 
ically down to zero temperature due to quantum fluc- 
tuations. The theory of quantum spin liquid is rather 
well developed, and it is expected that if such a state 
exists, it will have various exotic properties. For exam- 
ple, the low energy excitations may be objects which may 
carry spin ^ and no charge, called spinons. The spinons 
may be gapped or gapless and may obey either boson or 
fermion statistics. They will be accompanied by gauge 
fields, which may be of the U(l) or Z-i variety. In the past 
few years, several candidates for the quantum spin liq- 
uid have emerged. The best studied example is a family 
of organic compounds. The original k— (ET) 2 Cu2(CN) 3 
salt (abbreviated as ET) Q has recently been joined by 
a second material, [H the Pd(dmit)2(EtMe3Sb) which we 
shall refer to as dmit. Both materials are Mott insula- 
tors on an approximate triangular lattice with spin ^ per 
unit cell, but are not far from the Mott transition be- 
cause they become superconductor (ET) or metal (dmit) 
under modest pressure. There is no sign of magnetic or- 
dering down to 30 mK despite an exchange interaction 
J w 250 K. Both materials show a linear T coefficient of 
the specific heat at low temperatures, and a finite spin 
susceptibility. 01 The Wilson ratio is close to one, usually 
associated with metals and is highly unusual for an in- 
sulator. The thermal conductivity k is a good probe of 
these low lying excitations. Experiments on the ET salts 
indeed found a large contribution, but k/T is reduced 
below 0.3 K.[B[ On the other hand, recent experiments 
on dmit found that k/T extrapolates to a constant down 
to the lowest temperature. [1] These data strongly sup- 
port the picture that the low lying excitations are mobile 
fermionic particles, called spinons. 



Initial theoretical work pointed to a state where 
spinons form a Fermi surface and are coupled to U(l) 
gauge fields. Due to the proximity to the Mott tran- 
sition, ring exchange or charge fluctuations yield addi- 
tional terms to the Heisenberg model and help stabilize 
the spin liquid state. However, a peak in the specific 
heat around 6 K in ET and 4 K in dmit suggests a phase 
transition, which in the case of ET, has been confirmed 
by thermal expansion measurements. Q Furthermore, the 
nuclear spin relaxation rate 1/T±T shows a power law 
decrease below 1 K.[Io} These data led to the sugges- 
tion that the Fermi surface may be unstable to a pair- 
ing instability which nevertheless leaves a finite density 
of states intrinsically or due to impurities. ll| Thus the 
true ground state in the organic salts remains unknown 
at present. 

Two other examples, the Kagome compound 
ZnCu3(OH)gCl2 and the three dimensional hyper- 
Kagome Na^^Os also satisfy the condition of being 
spin liquids in that they do not show magnetic order and 
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both are characterized by gapless excitations. 
However, less detailed data are available and we shall 
focus our attention on the organics, even though the 
conceptual question we raise below will apply equally 
to these materials if fermionic spinons are found to 
be present. Spinons with Dirac sectra, however, may 
require a different treatment. 

In this paper we address two questions. First, how do 
the spinons couple to phonons, and, secondly, is there a 
way to unambiguously identify the pairing transition of 
spinons? As we shall see, the two questions are related 
because the attenuation of transverse sound turns out 
to be sensitive to the gauge magnetic field fluctuations 
and is a sensitive probe of the Meissner effect of gauge 
magnetic field at the onset of any pairing instability. 

The coupling of electrons to phonons is often discussed 
in terms of the screened Coulomb coupling between elec- 
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trons and nuclei and one may have the impression that 
the charge neutral spinon may couple differently. It turns 
out that in the long wavelength limit, the coupling ma- 
trix elements are exactly the same. This is because in 
this limit the coupling can be viewed as a distortion of 
the Fermi surface by the local stress of the unit cell and 
is the same whether the fermions are charged or not. Re- 
cently the spin phonon coupling was discused in terms of 
interactions mediated by gauge fields with the conclusion 
that the coupling is comparable to the electron phonon 
coupling. 1J] The present approach makes it clear that for 
a spinon Fermi surface, the result does not rely on gauge 
fields in the long wavelength limit. We next consider the 
relaxation of ultrasound due to spinons in parallel with 
the standard treatment of electron relaxation, including 
the onset of superconductivity. Differences and similari- 
ties will be pointed out. 

We begin with a derivation of the spinon phonon cou- 
pling following Blount's discussion of the electron phonon 
problem which was also used in Tsuneto's theory of ul- 
trasound attenuation. [ljj [JJI Blount's key insight is to 
transform to a frame moving with the lattice distortion. 
In a slight departure from Blount, we assume that the 
kinetic energy of the spinon is described by a mean field 
band Eo(k). Eo(k) can be a nearest neighbor tight bind- 
ing band, for instance. The unperturbed Hamiltonian 



H =E (p) + V imp (r'), 



(1) 



ijfp is the momentum operator, V[ mp de- 



where p - 

scribes disorder scattering which relaxes the spinon dis- 
tribution to the lattice, and r' refers to the laboratory 
frame. The sound wave is described by 5R(r',t) which 
is a smoothly varying function of r' such that it equals 
the displacement of the ions at the lattice points. The 
transformation to the moving frame r = r' — SR(r',i) 
is accomplished by a canonical transformation U = e lS 
where S = \{p- SR + SR-p). The transformed Hamilto- 



dU 

Ho + H 1 = UHU- 1 + i^-U- 1 
at 



(2) 



where Hq is the same as Eq.(l) with r' replaced by r, 
and keeping only first order in SR, 



H 1 =i[S,E (p)]+ l -Y t 



„ dS _„ 
at 



(3) 



Next we write SR ~ e l (i- r - ult ) where u) — v s q and v s 
is the sound velocity. We assume a slowly varying dis- 
placement, and compute the first term in Eq.([3]) to lowest 
order VSR and 2Jf. We find 



dSR . dSR 

H i = > j — — ppv a + ip ■ 



^-i dr a 

a/3 



at 



(4) 



where v a (p) — dEo/dp a is the electron velocity. Equa- 
tion (QJ is derived by formally expanding E(p) in power 
of p. The second term in Eq.@ is of order ujkp which is 
smaller than the first term by the ratio uikp / {qkpvp) — 
v s /vp and can be dropped. 

Now we introduce the phonons 



Hph = 



q,X.(7 



J q \a qX a q x 



(5) 



where A denotes the phonon branches with polarization 
i q \. Expanding SR in terms of the phonon coordinate 
and substituting in Q results in the spinon-phonon cou- 
pling term 



M kX (q)fl 



qa 



fka(a q x + a 1 _ qX ), (6) 



Mk\(q) = (k ■ £q\)[q ■ v(k)](2p ion uj q x 



,-1/2 



(7) 



where pi on is the ion mass density. The spinons are cou- 
pled to a gauge field a which initially has no dynamics be- 
cause it was introduced to enforce the the constraint of no 
double occupation. We shall approximate E(k) = k 2 /2m 
and the spinon kinetic energy is 



f 



(8) 



In addition we have the impurity scattering term H- lTCLV = 

Y,ia v ™p( r i)fl( r i)f<?( r i) and we consider non spin flip 
scattering only. We assume that impurity scattering 
gives an elastic scattering lifetime r and mean free path 
/ = vft for the spinons. 

From this point on we can discuss the sound attenua- 
tion in spin liquids in parallel with the theory for met- 
als and superconductors. Historically the first discussion 
was in the hydrodynamic regime valid for ql <C 1- [13, Ell 
The fermions are treated as a viscous medium subject to 
strain fields set up by the sound waves. This picture is 
clearly independent of the charge of the fermions and can 
directly carry over to the spinon case. Starting from the 
linearized Navier-Stokes equation, the sound wave relax- 
ation time t, is 
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PionVi 



T T = 

) ' a 
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PionVi 



TV 



(9) 



for the longitudinal and transverse sound, respectively, 
where rj, x are the shear and compressional viscosities. 
The sound attenuation constant a defined as the inverse 
of the phonon mean free path Z p h is the imaginary part 

of k = (oj/v s )(1 + iuJTs)^ 1 / 2 and given by a = in 
the limit ujt s <C I. The fermion viscosity is given by 
rj — j^N(0)m 2 v F T where N(0) is the density of states at 
the Fermi level and we can take x to be <C t]. 
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The hydrodynamic theory was extended by Pippard to 
all values of ql using a Boltzmann equation approach. [ijj 
Pippard pointed out that when ql > 1, the electrons de- 
velop local charge and current fluctuations for the longi- 
tudinal and transverse phonons which contribute signifi- 
cantly to the sound attenuation. Here we re-derive Pip- 
pard's results using a diagrammatic method, because it 
can readily be extended to the pairing case. Our method 
is simpler than the work of Tsuneto, who combined a dia- 
grammatic and Boltzmann approach. Since the diagram- 
matic treatment is not readily available in the literature, 
we provide the details in the supplementary material. 



Yl 0) = o2o/x * x\oo, = 



IX 1 ' = <"><>v x % x 




FIG. 1: Feynman digrams for phonon self energy. 



Let us first re-derive the results for metals. We com- 
pute the phonon self energy H(q, ui) due to the excita- 
tion of fermion particle-hole pairs. It is given by the 
diagrams shown in Fig.l. The bold solid line is the 
fermion Green's function with self energy due to im- 
purity scattering, G TCt{adv) (k,uj) = (u - £ fc ± i/lr)' 1 
where = k 2 /2m — [i. The electromagnetic (EM) 
field propagation is screened by repeated bubbles rep- 
resenting density or current fluctuations for longitudi- 
nal and transverse sound respectively. For longitudi- 
nal sound the Thomas Fermi screening length fc^p is 
much shorter than q~ 1 and it can be shown (see supp.) 
that the effect of screening is the same as calculating 
the unscreened bubble [Fig. 1(a)] with a coupling ma- 
trix Mk\ given by Eq.® with the trace subtracted, i.e., 
M kx (q) = [(k ■ e q \)(k ■ q) - \k 2 (q ■ i q \)}/m^2p ion u qX . 
This is the form suggested by Blount. [15] The standard 
results are reproduced, and the attenuation decreases in 
the superconducting state when the quasiparticles are 
gapped. In the rest of the paper we focus on the trans- 
verse sound. Here there is no vertex correction (impu- 
rity line across the bubbles) because the vertex given by 
Eq.flT]) is odd in the k component along q. describes 
the dissipation due to the creation of particle-hole exci- 
tation of fermions. In the particle-hole excitation 
creates current j which couples to the electromagnetic 
gauge field via j ■ A. The gauge field acquires a self en- 
ergy by coupling to current fluctuations, as shown in the 
third line in Fig.l. The resulting retarded photon prop- 



(10) 



agator (double wavy line) is given by 

£)EM _ 1 

rct iojcr±(q, ui) + ui 2 — c?q 2 

and c is the speed of light. For ultrasound lu 2 in Eq.(fT0|) 
can safely be ignored. Following Pippard's notation we 
write a±(q,uj) = qctq where cto = e 2 nr/m is the DC 
conductivity, and 

3 1 f 1 u n 

g = — Re[s 2 (a) - s (a)], s n (a) = — / du — — — , 
la 2i 7-i u + i/a 

where a = ql / (1 + ilot). We can safely assume lot <C 1 for 
the rest of the paper and set a — ql. Then g —> l-j^(ql) 2 
when gl C 1 and g —> ^-(ql)^ 1 when ql ^> 1. We see 
that in the clean limit (ql 3> 1) the first term in the 
denominator of Eq. (|10p is nothing but Landau damping 
N(0)uj/vFq while in the opposite limit (ql <C 1) it gives 
rise to the classical skin depth k^ 1 where k 2 — cjcto/c 2 . 

The diagrams are evaluated to give (see supplementary 
materials for details) 



LoN(0)k FSl (ql)-s 3 (ql) 



iql 



(11) 



On the other hand, n^' is proportional to the photon 
propagator, and depends on the relative size of the in- 
verse skin depth ho, q and I . Let us consider the 
case when c 2 q 2 <C Loa±(q,uj), in which case -Dj|t = 
(ioja±(q,Lu)) . This holds under the condition q <^ k if 
ql <C 1 and q 2 <C k^/(ql) if ql > 1. (We shall not discus 
the extreme clean case, q 2 ^> kg/ql while q <C fc , [l9[ 
because it is never attained for the spinon case.) We can 
show that Illg}. takes the remarkably simple form 



nS = — 9 -Imli°2 + 0(v s /v F ). 
9 

The ultrasound attenuation coefficient is given by 



Aim(nW + n«) 



nm 1 — g 

PionV s T g 



(12) 



(13) 



where the identity Si — s 3 = — y(l — g) + 0(v s /v F ) has 
been used, and n is the fermion density. Equations (fTB")) 
agree with Pippard's result derived using Boltzmann's 
equation. Using the limit g = 1 — j^(ql) 2 for ql <C 1, we 
can verify that the hydrodynamic limit is reproduced. 

Now consider the onset of superconductivity. II^ ) de- 
creases below T c due to the opening of the energy gap 
(see supp.), but nJ 1 ) is affected much more dramatically. 
Physically the onset of Meisner effect suppresses the mag- 
netic field fluctuations and II ^ ^ drops to zero rapidly be- 
low T c . Mathematically this is because a constant term 
e 2 n s (T)/m where n s is the superfluid density is added 
to the denominator of Eq. (fTUf and quickly dominates 
iwa±_(q,tS). Since II W is proportional to D™, it drops 
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rapidly to very small value. This is called the "rapid fall" 
and occurs over a millikelvin scale [20j j in clean samples. 
We note from Eq. (fT2|) that the fractional size of the drop 
is (1 — g) which is very small (~ j^(ql) 2 ) for <?Z <§; 1 but 
almost unity (1 — |^) in the clean limit of ql ^> 1. 

Next we turn our attention to the attenuation of trans- 
verse sound by spinons. The main difference is that the 
spinons and gauge fields are treated in 2D. Furthermore, 
the Maxwell term u 2 ~c 2 q 2 is missing in the photon prop- 
agator. The dynamics of the gauge field is generated by 
spinon current fluctuation and instead of Eq. (1101) we have 



Df et = 



1 



iua±(q,uj) - xi 2 



(14) 



where a± = gao, do = hit /to, n is the spinon density and 
X = 1/(247tto) is the Landau diamagnetism. 2jj Note 
that according to Eq.® the coupling constant to the 
gauge field has been set to unity instead of e. The factor 
g is calculated in 2D and is given by 



- 2 r , . , , . 1 

9 = - M a ) ~ t {a)\, t n (a) = — 
a li 



2 71 



d6 



cos 



(I 



cos 9 + i/a 



Once again, g can be considered a function of ql, g — 

(90 2 
4 



1 



for ql < 1 and g = ^ for ql > 1. 




FIG. 2: Schematic plot of the attenuation of transverse sound 
(for ql > 1) in a spin liquid normalized to qjv, the value 
in the Fermi liquid state. The "rapid fall" over a narrow 
temperature range AT at the transition to the spinon paired 
state is due to the Meissner effect of the gauge field. 

Just as in the EM case, we consider the case when 
Xq 2 <C uia±, i.e., Dj ct — (iwojj -1 and we conclude that 
a is given by Eg. (1131) with g replaced by g. In the Fermi 
liquid state, we can use this formula to get information 
on the spinon mass and lifetime r by studying the q (i.e., 
lu) dependence. At low temperature r is a constant dom- 
inated by disorder [(| while at finite temperature t may 
have interesting T dependence due to scattering by other 
spinons or phonons. If the spinons are paired, the rapid 
fall also occurs just below T c . Next we show that the con- 
dition XQ 2 ua± mentioned above is the only relevant 
limit. For ql > 1, u>aj_ is estimated to be (v s /vp)kp/m 
and will dominate XQ 2 as long as v s /vp 3> (q/kp) 2 . Since 



v s /vp 10 -3 , this condition is satisfied for most accessi- 
ble ultrasound frequencies. In the opposite case (ql <C 1) 
the condition is (v s /vp)ql >> {q/kp) 2 and is easier to 
violate. However, the rapid fall is very small in this case 
and difficult to detect and of little interest to us. 

Finally we can estimate the temperature range AT = 
T c — T of the rapid fall as sketched in Fig. 2. Since iwa± 
is replaced by iua — n s (T) /to in Eq. (|14[) below T c we 
find that Imni 1} = Imll^V (l + (n s (T)/mwcfog) 2 ) where 
IT 1 ^- are the values in the superconducting and normal 
states. For ql ^> 1, attenuation is dominated by IlW 
and we estimate AT as the temperature when ImFIs 1 "* 
has dropped half the normal state value. We assume the 
mean field (BCS) behavior n s {T) = 2n(AT/T c ) and we 
find 



AT 

~t7 



Vs 

vf' 



(15) 



We conclude that fermionic spinons in a fermion couple 
to phonons in a way which is identical to electron phonon 
coupling in the long wavelength limit. For ql ^ 1 the 
attenuation of transverse ultrasound is dominated by a 
component which is due to the fluctuations of transverse 
gauge fields. At the pairing transition of the spinon, this 
component is suppressed by the onset of the Meissner 
effect for the gauge field, leading to a rapid drop in at- 
tenuation in a very narrow temperature range below T c 
given by Eq. (fl"5|) . After the rapid fall, the attenuation 
is reduced in the usual way by the gapping of the quasi- 
particles. We believe the phenomenon of the rapid fall in 
the attenuation of transverse ultrasound gives a clear sig- 
nature of spinon pairing and the existence of U(l) gauge 
fields. 
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SUPPLEMENTARY MATERIAL 



Sound absorption and attenuation in a liquid 



Fluid viscosity will cause sound absorption and attenuation in a liquid. To discuss this effect, we begin with 
linearlized Navier-Stokes equation, 



du 



-Vp 



V + X ) V (V ' u ) — x V x u, 



(16) 



where p is the fluid density, u is the flow velocity, p is the pressure, rj and x are the shear and compressional viscosities 
respectivley. The instantaneous density p may be further written as p = po (1 + s), where s is a small fraction and po 
is a constant. Then equation of coninuity has the following form 



„ ds 



The acoustic pressure p is found to be of the form 



dp 
dp 



Pqs = Pov s s 



in terms of s and the sound velocity v s . So that we have a lossy wave equation 

5 \ 2 1 d 2 p 
1+Ts di) Vp= vlW> 



(17) 



(18) 



(19) 



with a relaxation time given by Eqs.(lO). If we assume monofrequency motion, the above wave quations is reduced 
to a lossy Helmholtz equation 



where k 



V 2 p + k 2 p = 0, 
The solution for a = — Imfc reads 



(20) 



1 UJ 

a = —= — 



1 + {ujT s y - 1 



1/2 



i + {uT S y 

In the limit lut s <C 1, we have the sound attenuation coefficient in Eq.(ll) 



(21) 



Screened spinon phonon coupling for longitudinal mode 

Without the loss of generality, the spinon phonon coupling Af^A (q) f° r longitudinal sound can be written in terms 
of spherical harmonic functions, 

MkA (q) = J/(fc, (?) + ai ™ Yl ™ ' ( 22 ) 

l>l,m 

,2 

where f(k,q) = . q , 9 is the angle between k and q, 6 is the azimuthal angle of k. For longitudinal mode, 

the Thomas Fermi screening length fc^p is much shorter than g -1 . We shall show that the monople part \f{k, q) will 
be screened by charge fluctuation, resulting in an effective coupling matrix AfkA (q, containing only higher order 
terms in multipole expansion (|22j) and a term of order of v s /vf- 

We firstly consider screening effect of the bare spinon density-density interaction L^ ' (q, u>), which gives rise to the 
screened interaction 



wmm - vwwwv 



M = 



^ 1} = 



■< 



FIG. 3: (a) Screening effect of spinon density-density interaction U, zigzag lines denote unscreened interaction U^ ' , double 
zigzag lines denote screened interaction U, the bubble is for the polarization function P. (b) The bare spinon phonon coupling, 
(c) Screened part of spinon phonon coupling. 



where P is the polarization function in an impure fermi system which is given by the bubble diagram in Fig. 3(a). 
The screening of spinon phonon coupling is very similar to that of U. These two can be evaluated through the same 
vertex function. We divide the screened spinon phonon coupling M^x (q, w) into two parts, 



M kA (q, w) = M kA (q) + (q, w) , 



where Afj^ in Matsubara frequency is given by 

M k\ (l) iuj n) = o 51 MkA ( q ) ^ ( k + ^' £ ( k > + ^™) ^ ( q > 7 (q: k , ik m ,iuj n + ik m ) 

as shown in Fig. 3(c). The vertex function 7 is the solution of the following equation, 
7(q,k, ik m ,iw n + ik ^ G (k'+q, ik m ) G (k', iu) n + ik m ) 



xT, 



k+q,k'+q 



(zfc m ) T k / k (zw„ + i 

h in 

)7 (q,k',ifc m , iuj n + ik 



where T k ' k is the transition matrix due to multiple impurity scattering and rii mp is the impurity density. 
To proceed we shall make use of the following approximations, 



£ a dv (k + q, e) G lc t (k, e + oj) 

/ 7: — : n F (e) Giet (k + q, e) <7ret (k, e + uj) 
J 2tti 



n F (e) £ adv (k + q, e) £ adv (k, e - u) 



qvp cos 9 + ui + i/r : 
dn F (£ k ) 

aCk ' 

dn F (£ k ) 



96 



For Z > 1, one sees that 



y dfiy; m (9, 4>) J 7£- n F (e) £rct(adv) (k + q, e) £ re t(adv) (k, e + w) = 0. 



With the help of optical sum rule and using the relation 

1 

27 ~~ ~ 7 

we find that 



'imp 



ImTl 



kk. 



71 (q, k, e, e + w) = 7 (q, k, e - 177, e + w + irj) = — rYT^ 

1 + s (a) / a 

12 (q, k, e, e + w) = 7(q,k,e + i?7,e + w + iry) = 1, 



(24) 



(25) 



(26) 

(27) 
(28) 
(29) 

(30) 



(31) 
(32) 
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where a — ql/(l + iu>r). Hereafter we shall assume ujt — ^ql <C 1 and set a = ql. So that the polarization function 



P can be evaluated as follows, 

P (q, iuj n ) = ^2^g (k + q, ik rn ) Q (k, iuj n + ik m ) 7 (q, k, ik m , iu n + ik m ) , 



k ik n 



Taking analytical continuation, we have 

/°° de 
—n F (6)72 (q,k,e,e + w) 
-00 2lT 



x [</ re t (k, e + u) Q let (k + q, e) - £ adv (k + q, e - u) G adv (k, e) 
de 



+2i 



-I 



2tt 



«f (e) 71 (q,k,e,e + w) 



x [-Q Tet (k, e + w) £ adv (k + q, e) + £ a dv (k + q, e - w) £ rot (k, e)] 
s (a) /a 



-JV(0) 



1 — iujt 



1 + s (a) I a 



When the Thomas Fermi screening length fc T p <g x , N (0) f/(°) (q,w) » 1, U (q, w) can be reduced to 



tf(q,w) 



JV(0) 



1 — iwr 



s a (a)/a 



l+s (a)/a 



Then we are ready to calculate (q), 



M kA (q> w ) = -s^2M kx (q)^2g(k + q,ik m )g(k,iu} n + ik m )jU(q,iuj n ) 



ik r) 



= \f (k, q) P (q, w) C/ (q, w) - 2iw/ (fc, g) £ (cos 2 - i 



/ ,oc de 


<9n^ (e) 


/-oo 2^ 


<9e 



C/ re t (k, e + w) t/adv (k + q, e) 71 U (q, w) 



/(*,«) 



1 .cjt S2 (a) — so (a) /3 

— [— ^ 

3 a 1 + (1 — iwr) sq (a) /a 



Therefore 



M k A (q) 



k 2 q 



ion^qA 



cos 6 — - + i 



1 , .D s s 2 (a) - s (a) /3 



3 vf 1 + (1 — iwr) sq (a) /a 



(33) 



(34) 



(35) 



(36) 



Thus we can calculate the sound attenuation constant through the unscreened bubble [Figl.(a)] with the traceless 
coupling matrix, 



MkA (q) = (k - q)(k ' g / A) ^ fc (q ' gqA) + O 



Vs_ 

vf 



(37) 



for longitudinal mode. 



Longitudinal sound attenuation 

Then we calculate longitudinal sound attenuation constant using the above spinon phonon coupling. In this case, 
we can express II (q, iu> n ) in terms of a vertex function T (q, k, ik m , iu> n + ik m ), 

II (q, iw n ) = - ^2 M kA (q) £ (k + q, ik m ) Q (k, iu n + ik m ) T (q, k, ik m , iu) n + ik m ) , (38) 



where the vertex function T obeys 

iuj n + ik m ) = M k \ (q) + ^ 5 (k'+q, «fc m ) £ (k', iw n + ik m ) 



^k+q,k'+ q {ikm) ?k'k (iu n + ik m ) T (q, k', ik m ,iuj n + ik m ) , 



and is shown in Fig. 4 diagramatically. 






+ o a 



FIG. 4: Vertex correction in longitudinal sound, the dash line across the vertex denotes impurity scattering. 

It is more convenient to write the analytical continuations of T as follows, 

Fi (q, k, e, e + ui) = T (q, k, e — iij, e + to + i-q) , 
r 2 (q, k, e, e + co) = T (q, k, e + i-q, e + u) + iff) . 

We find out the following two solutions, 



r x (q, k, e, e + u) = M kA (q) 



qk 2 F [s 2 (a) - |sp (a)] /a 



m v /2p 40 „iu q A l + s (a)/a 
r 2 (q,k,e,e + w) - M kA (q) . 
Thus the retarded phonon polarization function can be calculated as 

/°° de 
7T n F ( e ) I^ret (k, e + w) £ re t (k + q, e) T 2 (q, k, e, e + w) 

-^adv (k + q, e - w) CJadv (k, e) r 2 (q, k, e, e + w)] 

/°° de 
— n F (e)ri (q,k,e,e + w) 

x [-Giet (k, e + w) g adv (k + q, e) + £ a dv (k + q, e - w) £ rot (k, e)] . 
The imaginary part of the terms with T 2 and T 2 will vanish, resulting in 



r°° de 


drip (e) 




de 



ujN (0) fc| 
2mpi on v. 

Using the expressions of s n (o), 



F Re < s 4 (a) - -s 2 (a) + -s (a) 



Qret (k, e + w) C/adv (k + q, e) Ti (q, k, e, e + w) 

[s 2 Q) ~ | so (a)] 1 
a + sq (a) I 



s (a) = 
82(a) = -'- 

s 4 (a) = - 



—3a + a 3 + 3 tan 1 a 



3a 4 



we obtain the longitudinal sound attenuation constant, 

ujN(0)k% 1 



mpionvl 3ql 



q 2 l 2 tan" 1 (ql) 1 
ql — tan -1 (ql) 3 



q 2 l 2 tan" 1 (ql) 1 
ql — tan -1 (ql) 3 
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3 m 



is used. Eq. (|45j) coincides with Pippard's longitudinal sound attenuation 



where the density of spinons n 
constant. 

Similarly, we obtain the longitudinal sound attenuation constant in 2D by replacing s n (a) by t n (a) so on and so 
forth, 



^^Re <^ t 4 (a) - t 2 (a) + -t (a) 



nm 



[t 2 (a)-^o (a)]' 
a + t (a) 



q 2 l 2 



i - Vi + q 2 i 2 2 



(46) 



Transverse sound attenuation in an electron gas 

We proceed with transverse sound, where there is no vertex correction (impurity line across the bubbles) because 
the vertex is odd in the k component along q. The single bubble diagram n(°) can be evaluated as follows, 

/CO J 
— {[n F (e) - n F (e + u)] & et (k, e + u) G, dv (k + q, e) 

-n F (e) [0 ret (k, e + w) Q Tet (k + q, e) - £ adv (k + q, e - w) £ adv (k, e)]} . (47) 

Similarly to longitudinal mode, the imaginary part of GretQret and £J adv <?adv terms vanishes. So that 



/c 
-i 



^M4 Re[S2( 

wiV (0)k FSl (ql)-s 3 (ql) 



r°° de 


dnp (e) 


J-oo 2tt 




- S4 (a)} 





G ict (k, e + lo) G adv (k + q, e) 



Ap ion mv s 



iql 



where the relation 



s n (a) = - 



-s n -i (a) 



is used to derive the last line in Eq.(j48|) and we keep all the terms up to the leading order of v s /vf- 
Now we go on calculating IlW. Since the gauge field propagator is of the form 



H^ 1 ) (q, iu)„) can be written as, 
n«(q,iu; n ) 



fP-rri- 



D aj3 (q, iu) n ) = (5 a p - q a qp/q 2 ) D T (q_,iu n ) 



■ (q) Yl Q ( k + * lfc ™) Q ( k 



ct/3 k 



X] M k'A (q) (2fc Q + g Q ) (2fc^j + qp) (S a p - q a qp/q 2 ) D T (q,i< 

k' 

]T (k'+q, ^fc; n ) 5 (k', iw„ + ik'J 

, , # T (q, i0J n ) V M kA (q) A/k'A (q) [k • k' - (k • q) (k' • q)] 
p z m z A — s 



4e 2 



kk' 



^2 Q ( k + q> * fc m) 5 (k, iw„ + i/c m ) ^ (k'+q, ik' m ) Q (k', iw„ + ik' m ) , 



(48) 



(49) 



(50) 



(51) 



Summation over ik m and ik' leads to 



II« (q, iw„) = — D T (q, iw„) M kA (q) M k<A (q) [k • k' - (k • q) (k' • q)] 

kk' 

/OO 7^ 
— n F (e) [Q (k, e + iuj n ) A (k + q, e) + Q (k + q, e - iu>„) A (k, e)] 
-oo 7T 

r°° de' 

x / — n F (e') \Q (k', e' + iw n ) A (k'+q, e') + Q (k'+q, e' - iw„) A (k', e')] 

J-00 7T 

where the spectral function A (k, e) is given by 

G lct (k, e) - £ ad v (k, e) = 2ilm£ rct (k, e) = -iA (k, e) . 



Analytically continue iu) n —tw + iri, we have 

D 2 



nS (q, w) = (q, w) ]T M kA (q) M k , A (q) [k • k' - (k • q) (k' • q)] 



kk' 



f 00 de 

x / — n F (e) [£ ret (k, e + cj) A (k + q, e) + £ adv (k + q, e - w) A (k, e)] 

J-oo K 

r°° de' 

x / — n F (e') [& et (k' ; e' + w) A (k'+q, e') + g adv (k'+q, e' - ui) A (k', e')] 
After integrating over k and k', £ r ct<?rct and £ adv <? ad v terms vanish, we obtain 

(q, w) = -4^ r T ct (q, w) ]T (q) M k , A (q) [k • k' - (k • q) (k' • q)] 

kk' 

/OO J 
—n F (e) [-G iet (k, e + u) £ adv (k + q, e) + £ adv (k + q, e - w) C/ ret (k, e)] 
-oo * 
/•°° AJ 

x / — n F (e') [-Gret (k', e' + w) £ adv (k'+q, e') + S adv (k'+q, e' - w) £ re t (k', 

; 

-# r T ct (q, w) J] M kA (q) M k , A (q) [k • k' - (k • q) (k' • q)] 

»ret 

(k, e + oj) £ adv (k + q, e) 
G rct (k',e' + w)g adv (k'+q, e'). 



e 2 w 2 



kk' 



/■°° de 


<9n F (e) 


/-oo 


de 



r — 

J-oo K 



dn F (e') 



de? 



Therefore 



(q,a;) 



■•^ret (q^ W i ■ 



2m 2 
gv F 



4p lon mv s 



(o) 



- iF(q,w)ImII^(q,w) 



where 



F(q)W) = _e 2 4^(0) 



2m 2 gt> F 

For EM field the propagator Dj ct is given by Df™ (q, w) 



D Lt (q, [si (gZ) - s 3 (?0] • 
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Following Pippard's notation, we write a± (<?,w) = g&o, where oq — e 2 nr/m is the DC conductivity. To determine 
the factor g, we shall calculate cr± (q,u>) using the relation, 



o-± (q,u) = India (q,uj) , 



(58) 



where ILj± (q, iuj n ) = Jq d,Te lUnT Tlj± (q, r) and IIfj_ (q, t) = — (T r jj_ (q, r) • j j_ (— q, 0)) the retarded correlation func- 
tion. To be simple, we shall further count H± (q, oS) as the single bubble contribution. So that 



2e 2 1 2 

Hj_L (q, IWn) = o I k_ ( k • q) q] ^ ^ ( k + q> ^ ( k ; iuJ n + ik m ) ■ 

pm z 2 ^— ■ 



(59) 



Taking analytical continuation and then we have the imaginary part of (q, u>) 



Imn^(q,c) = -^Re5> 2 -(k.q) 2 ] f°° |l 



uje 2 k 2 F N(0)n ] 

27T77J 2 gup 



/■°° de 


drip (e) 






-s 2 (a)}. 





G Tet (k, e + w) C/ a dv (k + q, e) 



(60) 



Then we obtain 



a ± (q,u) 



e 2 k 2 F N (0) 
2m 2 qvp 
3cr 



Re [s (a) - s 2 (a)} 



2gZ 



S2 (ql) ~ s (ql)] , 



(61) 



where ctq = ^^m 2 *^''" * s use< ^- So that the factor g reads 



2ql 



[s 2 (ql) - s (ql)} 



(62) 



We find that g — > 1 — 2 ^ when gl « 1 and g — > when gZ ^> 1. We are interested in the case when c 2 q 2 <C 
wer_i_ (q,w), say, g <C g _1 / 2 fco, where fc^ 1 is the classical skin depth. This holds under the condition q <C ko if ql <C 1 
and g 2 <C fcg/ (?0 ^ ^ ^ ^ n this case ; 



D 



EM 



1 



icjcr^ (q, u>) ' 



(63) 



and 



Using the relation (|49|) . we have 



F(q,w) = - 



e 2 fcg^ N(0) gl [si (gi)-a3(g0] 
2m 2 g«F iwcrj_(g,w) 



^ [si (gZ) - s 3 (gZ)] 



(64) 



si (a) 
S3 (a) 



is (a) 



a 



i is 2 (a) 
3 a 



So that 



si (ql) - s 3 (ql) 



2i _ . s (ql) - s 2 (ql) 
" 3 1 ql 



~(1 -<)) + <>[- 
Vf 



01^ 



(65) 
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and 



So that wc obtain 



F = ^-^ + o(— ). (66) 



nS (q, «) = ^Imn^ (q, «) + O ) . ((»T , 



Thus the ultrasound attenuation cocffecicnt reads 



« = -2 Im( nS + nS) = -^^l. (68) 



TRANSVERSE SOUND ATTENUATION IN A SPIN LIQUID 

For a spin liquid, the spinons and gauge fields are treated in 2D. The gauge field propagator D T is replaced by 

Dj ct = , - 1 -, (69) 

where aj_ = g&o, <5o = nr/m, n is the spinon density and \ — l/(247rm) is the Landau diagramnetism. Note that the 
coupling constant to the gauge field has been set to unity instead of e. 
Similary, we have 

Imn j± (q, w) f K ' Re t (gZ) - i 2 gZ , (70) 

7TTTI qVF 



and 



crj.^w) = --Imft a (q,w) 

= f — —Re [t (ql)-t 2 (ql)} 

m z qvF 

2a 



So that the factor g reads 



ql 



We also find that g — > 1 — w hen gl<l and g ^ ^ when gZ >• 1. 
The imaginary part of fi^ is given by 



Re [t (g0 - *2 («01 ■ ( 71 ) 



^Re [t (ql) - t 2 («01 • ( 72 ) 



and is given by 



ImfiS(q> „ ) = ^ !l MzW), (73) 

2p ion mv s iql 



II«(q,u;) = (l + F)Imn^(q,w), (74) 
2i 

.9 



F = * [ti (gZ) - i 3 («/)] • (75) 



So that wc still can write as 



n« = l^Imn^ + ( ^). (76) 
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Onset of superconductivity 



For transeverse ultrasound, a rapid fall of the ratio a/ajq below T c due to IlW has been discussed in the paper in 
the clean limit ql ^ 1. We shall focus on II^ ^ term, which gives rise to the ratio or^/a^ and is given by 



-ImIT 



(0) 



(77) 



n^ ' will decrease below T c because of the opening of the energy gap. 

During the process of the attenuation of a sound wave in such a fermionic system, the sound wave passes its energy 
to fermions by fermion phonon coupling and fermions dissipate this extra energy due to disorder to complete the 
attenuation. Therefore the fermion life time r should be finite to relax the sound wave. On the other hand v s C«f, 
it is different from many other relaxation mechanisms. So that despite ql > 1, the physically relevant situation is 
reached by lot <C 1. It means we can not take r — > oo at first even in the clean limit because finite r is still seen for 
lot <C 1. Otherwise, we would obtain a II^ ) which does not depend on ql and conflicts with Eg. (14"51) . It is interesting 
that when screening is included by il^ -* + IT 1 ) , the result does not depend on r and appears to agree with the infinite 
r calculation (see Eq. trTTTO . It is only when screening is suppressed by Meissner effect that we can see this distinction. 
Below we shall discuss in the limit ojt <C 1 <C At, which is most interesting for us. Note that our results differ from 
those in early literatures (e.g., see J. R. Cullen and R. A. Ferrell, Phys. Rev. 146, 282 (1966).), since they worked in 
the clean limit by setting r — > oo which is not valid for ujt <C 1. 

To study the pairing state with finite r, we shall adopt the following impurity avaraged Green's function, 



Q (k, lLO n ) 



T (k, ILOr, 



(iQn) - £k _ A k (w„) 
A k (u) n ) 



(iu) n ) -£,1- A k (ui n ) 



A k (w„) = A k (l 



2r^ + A k ' 



2tV^+aT 



for superconducting states and replace the terms with QQ by QQ — J-J-. The results should coincide to Eq. (|4"5)l when 
A k — > 0. We do frequency summation firstly. 

S = -jj- Q (k, ik m + ibj n ) Q (k + q, ik m ) - T (k, ik m + iui n ) T (k + q, ik m ) 
f dz 

= 7T~ n F ( z ) $ ( k ' e + iuj ri) Q (k + q, z) - J"(k, e + iu) n ) J"(k + q, e)] 

Jc 27r * 

f°° de 

= 7T^ n F (e) {Q (k, e + iu n ) [G (k + q, e + irj) - Q (k + q, e - ir\)\ 

J-oo 2ttj 



+g (k + q, e - iuj n ) [Q (k, e + irj) - Q (k, e - irj)] 
-T (k, e + iio n ) [T (k + q, e + irj) - T (k + q, e - irj)] 
-T (k + q, e - iu n ) [F (k, e + irj) - T (k, e - ir/)]} . 



(78) 



The interval (—A, A), for which the Green's functions are continuous across the real e axis, can be excluded from the 
e integral. Taking iu> n — > to + irj, we have 



de 
2ni 



n F (e) {£ret (k, e + uj) [G let (k + q, e) - G adv (k + q, e)] 



-Sadv (k + q, e - w) [0 ret (k, e) - G adv (k, e)] 

- F ret (k, e + uj) [F let (k + q, e) - F adv (k + q, e)] 

-T adv (k + q, e - ui) [T rot (k, e) - T adv (k, e)]} , 



(79) 
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where 



Sret(adv) ( k > e ) 
?\ rct(adv) (k, e ) 
<?rct(adv) (k, e + Uj) 
•^"ret(adv) (k, £ + ^) 



t2 _ A 12 ' 
?k ^k± 



,'2 



_ rf _ A'2 ' 
k± ?k ^k± 



e k± ?k 



A" 2 ' 
^k± 



A" 
^k± 



112 
£ k± 



with 



e k± 



e 1± 



2r^ 2 -A 2 / 



AU = A k 1 ± 



2r V / i^Aj 



e k± 



= (e + w) 1± 



2Ty/(e + u) 2 -Al) ' 



A" 
A k± 



- A k 1 ± 



2rJ(6 + ^-A 2 , 



Then we perform integation over k in II' ) which can reduced to an integral over v = cos 6 and £ k - To do this, we 
write 5 ro t(adv) (k, e) and J" re t(adv) (k, e) in terms of the poles of £ k , 



^ret(adv) (k, e) 
-^rctfadv) (k, e) 
<?rct(adv) (k, 6 + Uj) 
•^ret(adv) ( k , e + Uj) 



4± + a 



c k± 



6 



c k± 



A' 2 
^k± 



2C k 



k± 



1 



+ 



1 



1 k± 



a: 



k± 



=k± 



k± 



A'2 
^k± 



//2 
£ k± 



_ f2 _ A »2 
Sk ^k± 



A^ ± 



c k± 



A "2 
^k± 



2Ck± 
2C k ' ± 
2Ck± 



c k± - a c k± + a 
i i 



Ck± - & Ck± + & 

1 1 



+ 



Ck± - & Ck± + & 
l l 



Ck± ^ ^ 



where 



: k± 



A' 2 
^k± 



C k± = 

c k ' ± = y^l - A k 2 ± = Ck' ± 
c k = ^ 2 - A l 



2t' 



We find that the following relations will be useful in dctcrming the imaginary part of , 



a: 



k± 



k± 



A k 



Ck 



•k± 



c k± 



AT 



HI 
<>k± 

A" 
^k± 

HI 
^ki 



€ + lu 
A k 



A^ 



Integrating with respect to £ k leads to 



/ 



— ^adv (k + q, e) Q Tet (k, e + w) 

d£k 4+q- + ^ + ?^ + £ k 



7T2 



2C k 



k+q 



2C k ' 



k+ 



Ck+q- 
1 



+ 



Ck - 9«F^ C k+ q- + Ck + O^F^ 



1 



( £ k+q- + 



k+q- 



: k+ 



C k ' 



k+; 



2C k+q _C k ' + c k '+-c k+q _ + ^ 
1 



+ 



HI 

c 'k+ 
1 



+ 



Ck ' Ck+ + Ck _ 

- ( e lc+q- - C k+q -) ( e k+ - Ck+) 



+ 



2Ck+ a _Ck4 



( £ k+q- C k 



k+q- 



E k+ 



+ c k '+) 



^k+q-'5k+ 
1 



HI _ H 

t 'k+ ^k 



k+q 



— qvpv 



k+q- 



Ck+q-j ( £ k+ Ck+) 



2C k + q -C k+ 



where 



C k+ + C k+q _ + ^ 
dCk 



/ 



2C k + q _C k + 



TF(CL ± ) 



k+q 



— qvpv 



etc. is used. Similarly, we have 



- J'adv (k + q, e) J" rct (k, e + w) 
A' A" 

^k+q-^k+ 



1 



A' A" 

^k+q-^k+ 



2C k+q _C k '+ C k '+ - C k+q _ + 2C k+q _C k '+ C k '+ - C k+q _ - ^ 

i A' A" i 



A' A" 

^ A k+q- /A k+ 



So that 



/ 



2C k+q _C k ' + C k '+ + C k+q _ + 2C k+q _C k '+ C k '+ + C k+q _ - OVf" ■ 

rfCk 



[Sadv (k + q, e) G let (k, e + w) - J" a dv (k + q, e) F rct (k, e + w)] 



1 ( e k+ q - + C k+q _ ) (e k+ + Ck+) 



A' A" 

ZA k+q- ZA k+ 



2C k+q _C k + 

1 



k+q- 



Ck+ - C k+q _ + 

Ck+q- ) ( £ k+ — Ck+) — A k+q _A^ + 



C" 

^ k 4 



c k 



k+q- 



qvpv 



2C k + q -C k+ 

i (4+ q --Ck+ q -)(^ + +c k ' + )-Aj c+q _Aj: 



k+ 



2Ck+q-Ck+ 



(4 



Ck+ "T S k + q - 



+ c k . 



qvpv 



k+q- 



k+q- 



l k+ 



CkV) 



A' A" 

ZA k+q- ZA k+ 



2Ck+q-Ck+ 



C k '+ + C k+q _ - qv F v 
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For the terms of (? ro t£ret — •^rct-^rct and G&dvG&dv — ^adv-^advj we need only to replace e'_ and by e' + and Q' + in the 
above, so on and so forth. Then we find that only the terms of (? a dv<?rct — -^adv-^rct contribute to the imaginary part 
T (0) 



of H Ic f Considering [M^\ (q)] is an even function of v, we have 

[M kA (q)] 2 [l + «(«+")- A k A ^ 



Imll£t = 2cj7riV(0)Re 



2w7riV (0) Re 



L />0O 

/ 

i J\& k \ 



de 



2ir 



de 



yj(e + uf - A 2 - ^ e 2 _ A 2 +q + i/T + qvpv 



dnp (e) 



de 



dv / — 



|A k | 



2tt 



9n F (e) 
de~ 



[M kA (q) 



1 + 



e -A k A 



^ A 2 ? - ^ - A k+q + i/r + qv F u 



For a s-wave pairing state, Irnll^ can be evaluated explicitly, 

Imn^ c 2c^iV(0)Rc jf — 



/■°° de 


dnp (e) 


/_' 


J A 2tt 


de 





^ [M kA (q)] 2 
_! qv F v + i/T 



The ratio a( ) / is given by 



= 2n F (A)ImII$ ) . 



„(0) 



a 



(0) 
N 



2n F (A) , 



(81) 



(82) 



(83) 



in the limit ujt <C 1 <C At. 



